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Interference in Lattice Networks 

Martin Haenggi 



Abstract 

Lattices are important as models for the node locations in wireless networks for two main reasons: (1) When 
network designers have control over the placement of the nodes, they often prefer a regular arrangement in a 
lattice for coverage and interference reasons. (2) If nodes are randomly distributed or mobile, good channel access 
schemes ensure that concurrent transmitters are regularly spaced, hence the locations of the transmitting nodes are 
well approximated by a lattice. In this paper, we introduce general interference bounding techniques that permit 
the derivation of tight closed-form upper and lower bounds for all lattice networks, and we present and analyze 
optimum or near-optimum channel access schemes for one-dimensional, square, and triangular lattices. 

■ Index Terms 

Wireless networks, interference, time-division multiple-access, geometry 

I. Introduction 



^ ^ " A. Motivation and contributions 

HH Wireless networks where nodes are arranged regularly in a lattice have advantages in terms of coverage, for 
^ ■ example in sensor networks or for cellular base stations, and in terms of interference, since it is much easier to 
i— 1| devise good channel access schemes than in networks where nodes are randomly deployed or mobile. Despite this 
advantage, relatively little work has focused on the interference characterization in such networks. Furthermore, 
^ ■ interference results for lattices also provide bounds on optimally scheduled general wireless networks, since the 
! goal of scheduling is to maximize the spacing between a receiver and its interfering transmitters, while maintaining 
^ ; a certain density of transmitters. 

Q ■ In this context, this paper makes three contributions: 

■^j- . • We introduce general bounding techniques for the interference in lattice networks. 

• We apply these bounds to transmitter-centric MACs (MAC schemes that schedule transmitters without con- 
sidering the location of their respective receivers, such as CSMA-type scheduling without RTS/CTS). In this 
case, transmitters form a lattice, but interference has to be measured at the receiver, where the interference is 
necessarily larger than if it were measured at the location of the desired transmitter. The transmitter-receiver 
distance r results in excess interference. We show that quadratic approximations of the form I(r) « J(o)+c ex r 2 , 
where J(o) is the interference at the desired transmitter and c ex r 2 is the excess interference due to the offset 
of the receiver, are highly accurate for small r. These results are relevant for well scheduled wireless networks 
with arbitrary node distribution. 

• For networks where all nodes form a lattice, we analyze and compare the interference of different TDMA 
schedulers, and we provide schemes that are very close to optimum (if not optimum) for one-dimensional, 
square and triangular networks. In the one-dimensional case, we also provide results on the achievable rate 
and the transport capacity. 
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B. Related work 

While a growing body of work studies interference in random networks (see, e.g., [1], [2] and references therein), 
only few papers have addressed the issue of interference in lattice networks. In [3], bounds on the interference 
in triangular networks were derived using a relatively crude upper bound on the Riemann zeta function that is 
within 25% of the true value for the range of 2 to 4. We will derive a much tighter bound that is within 1.3%. 
A TDMA scheduling scheme for square lattices that is optimum for the case where the density of concurrent 
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transmitters is 1/4 is suggested in [4]. Here, we provide near-optimum scheduling schemes for any density 1/m 2 , 
m G N. The interference distribution in one-dimensional networks with Rayleigh fading is analyzed in [5] for the 
case where all nodes transmit, and [6] derives outage results an throughput-optimum TDMA schedulers for the 
same type of network. Finally, the single-hop throughput for two-dimensional lattice networks with Rayleigh fading 
is approximated in [7]. For non-fading channels, [8] provides throughput results for general TDMA schemes in 
two-dimensional lattice networks. The interference is expressed using complicated infinite double sums (that are 
evaluated numerically), for which we will present tight bounds. 

The network models in this paper are entirely deterministic, although, of course, the interference results derived 
correspond to the expected interference in fading channels. 

The paper is organized as follows. We first introduce the bounding techniques (Section II), followed by three 
application sections that discuss one-dimensional networks (Section III), square lattice networks (Section IV), and 
triangular networks (Section V). Conclusions are drawn in Section IV. 

II. Bounding Techniques 
In this section we introduce the basic techniques that will be used to bound the interference. 



A. Upper bounds 

Let L C R d be a <i-dimensional lattice, i.e., 



L^{x = Gu: ueZ d }, 



(1) 



where G G R dxd is the generator matrix. It is assumed that detG ^ to exclude degenerate cases. Important 
cases include the square integer and the triangular lattice in two dimensions, both with nearest-neighbor distance 
1: 
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(2) 



The lattice L has the properties that it includes the origin o = (0, . . . , 0) and that each lattice point is centered in 
its Voronoi cell, i.e., if VL(x) is the Voronoi cell of lattice point x G L and U x is a uniformly randomly distributed 
random variable on V^(x), then EU X = x. The volume of each Voronoi cell is V = | detG|, and the density of 
the lattice is A = V^ 1 (points per unit volume). 

Let L !o = L \ {o} 1 . Then the interference at point z is defined as 

I(z) — £{x — z) , 



where £(x) : R d -> R 
a function £'(r): R + 



is the path loss function, assumed monotonically decreasing and isotropic, i.e., there exists 



such that £' 



= £(x), where || • || is the standard Euclidean norm. Further assume 



that £'{r) = o(r~ d ) as r — > oo (otherwise the interference is infinite) and that £'{r) has a convex tail, i.e., there is 
a finite radius, defined as 

r c = inf{r > | £'(p) convex for p ^ r} . 

All commonly and rarely used path loss functions satisfy these properties. The origin is assumed to be the desired 
transmitter, so it does not contribute to the interference. The distance \\z\\ is restricted to values that ensure that z 
is not too close to an interferer, i.e., \\z\\ < \\gi\\, for 1 < i < d, where g, L are the column vectors that constitute G. 

The first theorem exploits the tail convexity of the path loss function to yield an upper bound on the interference. 
Let b x (r) be the d-dimensional ball of radius r centered at x. 

Theorem 1 (Voronoi upper bound) Without offset: Let L be a subset o/L !o such that there is no overlap between 
any of the Voronoi cells of the points in L and b {r c ). Then 



lip) < £ £{x) + i£ / 



£(y)dy. 



(3) 



xeV°\L 



eL JVL(x) 



'The superscript ! ° is borrowed from stochastic geometry, where it denotes a random point set, conditioned on having a point at the origin 
but excluding that point. 



3 



With offset: Let L be a subset of L !o such that L does not include any points whose Voronoi cells overlap with 
b (r c + IMD- Then 

xeV°\L xeL J W 

Proof: We only need to prove (4). Since £'(r) is convex for r > r c and the Voronoi cells of all points in L 
lie outside b Q (r c + ||z||), it follows from Jensen's inequality that £(x — z) = £(E(U X ) — z) ^ E(£(U X — z)) for all 
x G L, where U x is uniformly distributed on V^(x). ■ 
Next we state a corollary that applies to two-dimensional lattices and provides a bound that is simple to evaluate, 
as it is based on a radial outer bound on the integration region. 

Corollary 2 (Radial bound for two-dimensional lattices) Let p x = mm yeVL ^ \\y\\ be the (minimum) distance 
of the Voronoi cell of point x to the origin, and define 7Z Z = {p x : x G L !o , p x ^ r c + \\z\\}. Without offset: For 
any 6 1Z , 

I(o) < e ( x ) + 17 / r ^'( r ) dr • ( 5 ) 

_w .„ " Jr h 



x€L u 
p^<r b 



With offset: For any r^^Tl 



I(z) < Y £{x-z) + — r£'(r)dr . 

^ i„ ' Jr b -\\z\\ 



(6) 



Px<Tb 



Proof: No offset: Since the union of all Voronoi cells of the points not included in the sum is a strict subset 
of R 2 \ b (ri,), (5) follows from Theorem 1. With offset: The change from o to z means that £(x) is to be replaced 
by £{x — z) in the sum and the integral. The integral in (5) is taken over M 2 \ b (ri,), thus, with the offset, the 
domain of integration is E 2 \ b- z (ri,). Since this integral may be tricky to calculate, we replace it conservatively 
with M 2 \ b (rb — \\z\\) for the upper bound. ■ 
Remarks: 

1) For the case where the path loss law £(x) has a power law tail with exponent a, the integral evaluates to 
27r6 2 ~ a /(a — 2), where b is the lower integration bound. 

2) The reason why r& is restricted to the set 1Z is that if is not the smallest distance of a Voronoi cell, then r& 
can always be increased to the next larger smallest distance without changing the set over which the sum is 
taken, thereby yielding a better bound. In other words, for any set L chosen in Theorem 1, the corresponding 
rf, here should be the largest radius such that b a (rb) does not overlap with any of the Voronoi cells of the 
points in L. 

3) The integral in (5) can be interpreted as the mean interference stemming from a Poisson point process of 
intensity A(r) = l{ r>ro } /V. As one would expect, outside of a certain radius from the origin, only the intensity 
of the point process matters, rather than its higher-order statistics (such as second-moment measures). The 
corollary provides a bound on this radius such that the interference in the lattice is upper bounded by the 
interference in the Poisson point process of the same intensity. 



B. Lower bounds 

To obtain lower bounds, the distances of two or more points (in the convex region of the path loss function) are 
replaced by their averages. Evaluating the path loss at this average point, multiplied by the number of points that 
the average was taken over, yields a lower bound. This technique is not restricted to lattices. 

Theorem 3 (Lower bound for general point sets) Let L C K rf be an arbitrary discrete set of points that includes 
the origin, and let L'° = L \ {o}. For any LcLn b c Q (r c + \\z\\), let xl = (1/\L\) YlxeL x ^ e avera S e °f tne 
points in L. Then 

I(z)> Y, Kx-z) + \L\£{x L -z). (7) 

x£V°\L 
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Proof: This follows again from Jensen's inequality. Moving any set of lattice points (outside distance r c from 
the origin) to their average yields a lower bound due to the convexity of £ in this regime. ■ 
In general, this theorem will be applied repeatedly; for example, many pairs of points are formed and replaced by 
their averages. Better bounds can be expected if \L\ is kept small, and if the points in L are located nearby. 



C. Bounds on zeta function 

For power law path-loss functions, the interference in lattice networks can often be expressed using the standard 
zeta function (without offset) or the generalized zeta function (with offset). The following lemma provides tight 
yet simple closed-form bounds on the standard and generalized Riemann zeta function. 

Lemma 1 The generalized Riemann zeta function 2 

oo 

C(a,l-z) 4 a >l, \z\<l, 

k=l 

is tightly upper bounded by 

C(a, 1 - z) < (1 - z)' a + ^ >— . (8) 

a — 1 

and lower bounded by 

C(a, 1 - z) > (1 - z)- a + ^ 1— . (9) 

a — I 

For the standard zeta function Q(a) = ((a, 1), an alternative upper bound that is even tighter for a > 2 is 

C(«) £ " 7 m _ a - (10) 

a — l — (a — 1)2 a 

and a good lower bound is 

C(a) ^ 6a _ 3 a_ 2 a_ l - (H) 

Proof: The first bound (8) is an application of Theorem 1 with L = N and L = N\ {1}. Let be a random 
variable that is uniformly randomly distributed in [k — 1/2 — z, k + 1/2 — z). From Jensen's inequality, it follows 
that (EX k )~ a < K(X k a ) for all k > 1/2 + z and a > 0. Expressing the zeta function as 

n oo 

C(a,l-z) = J>-z)- Q + E ( E ^)~° 

fe=l fc=n+l 

for arbitrary n G N and upper bounding each term (EJi) _a by E(X^ a ), we obtain the bounds 



n POO 

((a, 1 - z) < ^(k - z)~ a + / x- a dx , 



iiGN, 

fc=1 ^n+l/2- 2 

which are increasingly tight as n grows. (8) is the bound for n = 1. For the lower bound (9), we change the support 
of all Xfc to [A; — z, k + 1 — z\. In this case, since X^ ^ k — z, we have (At — > E(X A T a ) and thus 



7! 



DC 



C(a,l-z) >^(k-z)~ a + / a;- a dx, n G N . 

The bound (9) is the bound for n = 1. Regarding (8), we consider the special case z = 1/2. It is straightforward 
to show that 

C(a,l/2) = (2 a -l)C(a), (12) 

2 also called Hurwitz zeta function 
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which, together with (8) for z = 1/2 yields (2 a - l)C(a) S2 a + l/(a - 1), which proves (10). The lower bound 
(11) stems from an application of Theorem 3 (with z = 0) and 

L = \J {-6k - 1, -6k + 3,6k- 3, 6k - 1} (13) 

fceN 

and averaging on the distances of the points — 6k — 1 and 6k — 1, /.<?., replacing them by two points at distance 
6k. Hence we compare the two sums 

S! = ^(3& - l {fceven} r a + (3* + l {fce ven } r a ! S 2 = 2^(3k)~ a , 

fceN fceN 
where S2 is obtained from Si by averaging over the two arguments 3k — 1 and 3k + 1 whenever fc is even. The 
terms for odd k are the same in S\ and S2. From Theorem 3, we have S\ > S2, where 

s 1 = c(a) - 1 + j>£r a - E( 2A; )" a - Ew = + - 2- a - 6- Q ) - 1 , 

ken fceN fceN 

and 5 2 = 2 • 3~ Q C(a)- ■ 

Remarks: 

1) For a < 2 (and z = 0), the first bound (8) is tighter than the second one (10). They are equal at a = 2, 
which is the value of a where the first bound is loosest. The difference to the actual value ((2) = tt 2 /6 is 
(10 — 7r 2 )/6 « 0.0217. The two bounds (10) and (11) are so tight that when plotting both curves over a range 
[cti, a u ] with a u — a\ > 1 they appear as one. For a I 1, the upper bound is tighter, while for practical values 
of a, the lower bound is. 

2) A simpler and looser lower bound is 

C(«) > 



2 a - 2 ' 

obtained from summing over 1, 2, 4, 4, 6, 6, 8, 8, . . . instead of N. 

III. One-dimensional Lattice Networks 

A. Interference with offset 

Here we assume that the desired transmitter is located at the origin o, the receiver at position < z < 1, and 
interferers at positions Z \ {o}. Unless otherwise noted, we use the standard power path loss law l(x) = ||x|| _a . 
Since I(z) = ((a, 1 — z) + ((a, 1 + z), we can apply Lemma 1 to obtain a tight closed-form upper bound: 



I(z) < (1 - z)~ a + (1 + z)~ a + ' 



a 



2- Z ) + {-2 +Z 



(14) 



Due to the symmetry of the arrangement, I(z) is an even function. For integer a, the bound is a rational function in 
z 2 , with the numerator polynomial of degree 3a — 1 — l{ a even }, and the denominator polynomial of degree 4a — 2. 
A quadratic Taylor expansion at z = gives a good approximation for small z (in this case, z < 1/4 is small 
enough for a very good approximation): 

I(z) w 2C(a) + c cx z 2 , 

where c cx = I"(0) = (l/2)d 2 I(z)/(dz) 2 \ z= o is the (transmitter-receiver) offset coefficient. We obtain 

a+l\ 



c cx = a 2 + a ^1 + (j: J j . (15) 

For all practical values of a, a good lower bound is c cx > a 2 + a, and a good approximation is c cx w a 2 + a + 1/2. 
It is also possible to give a simple yet accurate approximation for the fourth-order coefficient: 

4x } ~ + \a 3 + a 2 + l -a 

This approximation for cix is within 1% for all a ^ 1. For practical a and small z, the quadratic term always 
dominates the fourth-order one. 
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Eqn. (14) can be used to find the displacement z that maximizes the transport capacity T(z) = zlog 2 (l + 
z~ a /I(z)), where z is the link distance, and the logarithmic term is the rate R(z). Since 1 + z~ a / 'I(z) is a rational 
function for integer a, a series expansion of the form R(z) « c\ + c 2 log z + c^z 2 can easily be found. Multiplying 
by z and setting the derivative to zero leads to solving an equation of the form 1 — 2 log z/a + bz 2 = 0, which 
yields z opt = exp(a/2 — W(a6e a ))/2), where W is the Lambert W function. For a = 2, this yields z op t = 0.223, 
while a numerical investigation yields z opt = 0.224. For a = 4, the gap is slightly larger: The analytical result is 
z pt = 0.217, and the numerical one is z opt = 0.222. For non-integer a, similar series expansions can be derived 
by substituting x = z a and optimizing x. It can be concluded that the optimum value of z is between 0.2 and 0.25 
and rather insensitive to the value of a. 



B. Application to TDMA scheduling patterns 

1 ) Separating transmitters: Let the nodes at positions mZ be transmitters, each one transmitting to the next node 
on the right. Focusing on the desired transmitter at the origin and its receiver at location 1, the other transmitters 
at positions K. = mZ \ {o} are interferers. This is an m-phase TDMA pattern, since it will require m time slots to 
give each node a transmission opportunity. The interference at location 1 can be expressed as 

I m = ^2\k-l\- a =m- a I(l/m), m>l. 



With Lemma 1 and 



A (m-l) Q + (m + l) a {bm + I) 01 " 1 + {bm - I)*- 1 
G(m, b) — —75 — h 



we have 

The maximum achievable rate is 



(m 2 — l) a m(a — l)(b 2 m 2 

G(m,2) < I m < G(m, 3/2). 



(m a+2 \ 
1+ 2C(a)„,> + eJ - <16) 

where c cx is given by (15), and the throughput is T m = R m /m. The derivative dT m /dm indicates that the optimum 
to is always to = 4 or to = 5, with T4 T5 for all a 6 (1, 20]. Due to the smaller power consumption, m = 5 is 
preferred in practice even if m = 4 yielded a slightly higher throughput. So, interestingly, irrespective of the path 
loss exponent, m = 5 is an excellent choice — assuming the rate of transmission is adjusted to the interference 
present for that path loss exponent. A first-order Taylor expansion at a = 2 yields T5 w 0.60 + 0.53(a — 2), which 
is a very good approximation for the practical range of a. 

The throughput R m /m is identical to the transport capacity zR(z) considered in the previous subsection when 
z = 1/m. Thus the result that m = 4 or m = 5 are optimum for all a also follows from the fact that z opt is 
between 1/4 and 1/5 for all a. 

The optimum m for one-dimensional networks with Rayleigh fading was derived in [6]. It was found that 
"iopt = 3 for a > 2. So in the fading case, transmitters can be packed more closely. 

This m-phase TDMA scheme suffers from an imbalance in the interference from the left and right; nonetheless, 
it is optimum for unidirectional traffic. For bidirectional traffic, a better TDMA scheme exists, as discussed next. 

2) Balanced TDMA schemes for bidirectional traffic: Here we analyze TDMA schemes for bidirectional traffic, 
which achieve higher throughput than in the unidirectional case, since it is possible to balance interference from the 
two nearest interferers. These MAC schemes separate the receivers from their interfering transmitter in an optimum 
way. In other words, instead of placing transmitters optimally, these schemes place transmit-receiver pairs optimally. 

We use the following notation to define transmission patterns: Transmitters are indicated by T, with an arrow 
on top to indicate on which side their receiver sits, and receivers are denoted by R. As a shortcut, R fe indicates a 
sequence of k receiving nodes. The period of the transmission scheme is indicated by a subscript. For example, the 
unidirectional TDMA scheme with m = 3 is described as (T*RR)3 or (T*R 2 )3, and for general m, (T^R m_1 ) m . 

In the balanced schemes, transmitters with their receivers to the right alternate with transmitters with their 
receivers to the left. For m = 3, for example, the pattern is (T^RRRT^Rje. This has a period of 6, but since two 
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m=2 X — *0 CK — X X — >Q CX — X 



any m X-^O ... O-X O ... O X-^O ... O— 
> „ ' > „ ' » v ; 

m m-2 m 

Fig. 1. Balanced TDMA schemes for m = 2 and general m. The two nearest interferers are at distances m (to the left) and in (to the 
right). 



nodes transmit in a group of 6 nodes, it take 3 time slots for each node to transmit once. For general m, the scheme 
is 2 m - These balanced transmission patterns are illustrated in Fig. 1. 

Let m = 2 and the receiver under consideration be located at the origin and listening to the transmitter at — 1, 
i.e., we focus on the underlined receiver in the pattern The interferers are located at 

/C 2 = (J {— 4A — 1, — 4Ar + 2,4fe — 1, 4A; — 2} , 
ken 

and the interference is 

h = 1*1"° = ^ ~ 1 + X> fc )~ a " 2 E W = C(«)(l + 2" a - 2 • 4- a ) - 1 . 

keK fcgN fcGN 

^2 can be bounded by bounding the zeta function using Lemma 1. Alternatively, a lower bound is obtained by 
replacing the interferers at positions 4k — 1 and —4k — 1 by two at position 4k and applying Theorem 3: 



2 • 3° 



1 



For m = 3, 



/C 3 = |J {-6fc - 1, -6k + 3,6k- 3, 6fc - 1} . 

fcGN 

The exact expression is 

i 3 = J2 \ k \~ a = c(«) - 1 + I]( 3fc r Q - (2*o~ a + ( 6fc r a = c(«)(i + - T a - 6~ a ) - 1 

Using the lower bound for £ from Lemma 1, we obtain 



>a+l 



/ 3 > 



6 a - 3° - 2 a - 1 

Generalizing to arbitrary m > 1, there are interferers at /C m = {/c G N | —2km — I, —2km + m, 2km — m, 2km — 1}. 
Replacing each pair —2km — 1 and 2km — 1 by two interferers at 2km, we have 

7 m >2^(mfe)- Q = 2m~ a ((a), 
ken 

which can be (further) lower bounded by any of the lower bounds on the zeta function from Lemma 1, e.g., 



I m > 



2 • 6 C 



> 



1 + 



1 



m a (6 a - 3° - 2 a - 1) m a [ 2 a - 2_ 
To find a strict upper bound on I m , we replace the elements 2km — 1 in )C m by 2km — m + 1, which yields 

(a) 



ken 

= 2(2m)- a (k - l/2)- a + m- a (k + l/m)- a 



ken 



= 2{2m)- a Q{a, 1/2) + m~ a Q{a, -l/m) 



(b) 1 , 
< (2 + 

m°' 



)1— a 



a - 1 



+ (m + l)- Q + 



m(a — 1) \.3m + 2 



a-l 
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a=5 






a=4 






a=3 






a=2 




>- — — 



' 5 2 2.5 3 3.5 4 4.5 5 5.5 6 



m 

Fig. 2. Throughput as a function of path loss exponent a and TDMA parameter m for balanced schemes in one-dimensional networks. 
The dashed lines are the approximations log 2 (l + m a C(a))/m. 



where (a) is strict for m > 2 and (6) follows from (8) in Lemma 1. Since the lower bound is tighter and reveals 
that I m is essentially proportional to mr a , we proceed with the lower bound to find the throughput-optimum m. 
The SIR is I m l = m a C{a), where, from (11), 

C(a) = ^ « \ (1 - 2- - 3- - 6"°) . (17) 

With R m = log 2 (l + m a C{a)), the throughput T m = R m /m is maximized at 

_ W(a) + a 
opt C(a)W(a) ' 

(relaxing the integer constraint on m), where a = — ae _Q and W is the Lambert W function. The (real) m opt G [3, 4] 
for a < 7. Comparing T3 and T4 shows very little difference, and again power consumption favors T4. A second- 
order expansion of T4 as a function of a is T^a) « 0.64 + 0.59(a — 2). This is better than in the unidirectional 
case, but the difference is not large. Noting that 4 a C(a) w (4 Q — 2 Q )/2 for practical a, the optimum rate at m = 4 
is i?4 a — 1 + log 2 (2 a — 1), or just a bit less than 2a — 1. Fig. 2 shows the throughput as a function of m and 
a, together with the approximation log 2 (l + m a C(a))/m, where (7(a) is given in (17). 

IV. Square Lattices 

In this section, we consider square lattices L = 1? . In the case of the power law path loss, we may use results 
on lattice sums, see, e.g., [9], to express the interference at the origin: 

I{o) = £ \\x\\~ a = 4C(a/2)0(a/2) , (18) 

where 

«*>*£^ 

is the Dirichlet beta function. /3(2) is Catalan's constant K = 0.916. So, for a = 4, / = 27r 2 i-C/3 « 6.03. These 
expressions are not closed-form and are restricted to the origin, and thus do not provide much insight into the 
behavior of the interference. We apply the techniques described in Section II to derive tight closed-form upper and 
lower bounds. The path loss law assumed is again £(x) = ||x||~ a . 
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A. Lower interference bounds 

To obtain a good approximation on (18), we group the nodes into (square) "rings" of increasing distances. In 
the k-th ring, there are 8A; nodes at distances between k and \[2k. Let be the average of the distances of the 
k-th ring of nodes. Then, taking the sum over the nearest ring separately, 



I(o) >4(l + 2^/ 2 ) + 8^D A 



k=2 



by Theorem 3. Replacing by an upper bound ck > still provides a lower bound. A simple bound is 
D k < k(l + V2)/2. This is an upper bound since 

fc(l + vg) ^ 1 ^j k2 + . 2 + ^ fc2 - (fc _ ^ y ^ if ^ k _ 

So choosing c = (1 + \/2)/2 ps 1.207 yields the bound 

7(o)>4(l + 2-^) + 8 - 2 7« a -^- 1 > . 

(1 + V2) a 

For a = 4, this is about 5.76. The exact value is 6.03. Using the simple lower bound ((a — 1) — 1 > 2/(2° — 4), 
a > 2, 

o«+4 

1(6) > 4(1 + 2~ Q / 2 ) + — . 

V ' V ' (2«-4)(l + V2) Q 

The coefficient c can be sharpened. The average distance is given by 

fc-i 



D k = — 

k 8k 



4k + AkV2 + 8k^ \/l + J 2 A 2 



Noting that 



we obtain the bound 



For k > 1, Dfc ^ ck for 



V Vl + i 2 A 2 < (k - 1) / \/l + z 2 cb , 

D fc < -(l + V^+(fc-l)(V2-log(v^-l))) • 
c4^ + i(l-log(^-l))« 1.1775, 



which yields the sharper lower bound 

oo 

J(o) > 4(1 + 2~ a / 2 ) + 8k(ck)- a = 4(1 + 2~ Q / 2 ) + 8c~ a (((a - 1) - 1) 

k=2 

Lower bounding the zeta function using (11) in Lemma 1 yields the closed-form bound 

1(6) > 4(1 + 2~ a / 2 ) + 8c 

For a = 4, this is 5.84. 



t 3"- 1 + 2"- 1 + 1 

ga-l _ 30-1 _ 2<*-l - 1 ' 



(19) 
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B. Upper interference bounds 

Let Xij be a uniformly randomly distributed variable on [i — 1/2, i + 1/2] x [j — 1/2, j + 1/2]. Then for any 
Z C L !o , the lattice sum (18) is bounded by 

i(o)= E (^-r Q ^ E E (^-r Q + E E K7) 



E E (^-r Q +^ikir Q dx, 



where Z c = L !o \ Z and 



A = M 2 \ J it 1 " V2,* + 1/2] x [j - 1/2, j + 1/2]} 
\CO')ez 

where © is the morphological dilation operator and the structuring element is the unit square centered at o, i.e., 
un = Vl(o) = {x € M 2 : \\x\\oo ^ 1/2}. The bound is tighter if Z includes the points closer to the origin. For 
example Z = {x G I? : < ||x|| ^ \/2}. This gives the bound 



I(o) < 4(1 + 2~ a/2 ) + f \\x\\- a dx. 

Jm. 2 \\0.3/2] 2 



'R 2 \[0,3/2] 2 

This integral does not have a solution for general a. For rational a, it can be expressed using hypergeometric 
functions, and for integer a, it simplifies to 

a = 3: 8^2/3; a = 4: 2(2 + ^)/9; a = 5 : 80^2/243. 

So for a = 4, we get I(o) < 6.14, which is within 1.5% of the exact value. 

Changing to polar coordinates yields simpler closed-form bounds. As an application of Corollary 2, we outer- 
bound the area A by M 2 \ b Q (b) with 6 = 3/2. This gives 

I(o) < 4(1 + 2-/ 2 ) + 27r(3/2) o 2 "" , (20) 

a — 2 

which is a decent bound for a > 3. For a = 4, 1(6) < 6.40. Better bounds can be obtained by using a larger 
radius in Cor. 2. Adding 12 more nodes in the direct calculation and increasing = 3/y/2, we obtain 

1(6) < 4(1 + 2~ Q / 2 + 2' a + 2 • 5 -«/2) Mg/v^! . (2i) 

a — 2 

This situation is illustrated in Fig. 3, where the circle has radius r&. The lower bound (19) and the tighter upper 
bound (21) are shown in Fig. 4. 

C. Transmitter-receiver offset 

Let I(z) be the interference measured at location z, with \\z\\ < 1. First we consider the case where z lies on 
the positive real axis, i.e., z = (\\z\\, 0). Let r = ||z||. Let Is denote the interference from the 8 nearest interferers. 
From elementary geometry, the distance x to interferer (1,1) is given by x 2 = 2 — r(2 — r), and the distance x' to 
(—1, —1) is given by x' 2 = 2 + r(2 + r). Hence 

I 8 (r) = (1 - r)~ a + (1 + r)- a + 2(2 - r(2 - r))~ a ' 2 + 2(2 + r(2 + r))~ a / 2 + 2(1 + r 2 )" a / 2 
= 4(1 + 2~ a ' 2 ) + q 2 (1 + 2~ Q / 2 ~ 1 ) r 2 + 0(r 4 ) . 

\ / V J 

7 8 (o) c s 

A simple lower bound for 2 < a < 12 on the excess interference coefficient is cs > a 2 + 1. Including the 
interference from all other nodes: 



I(z)<I 8 (r)+ [ 



\\x — z\\ a dx 

K 2 + \[0,3/2] 2 
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Fig. 3. Voronoi cells for square lattice. The circle of radius 3/v2 indicates the radial bound for the integration for the bound (21). The 
interference from the 20 nearest points, marked by a cross, is summed up directly, while the interference from the other nodes is upper 
bounded by integrating H^H -0 over the outside of the circle. 

1 8 1 1 1 1 1 1 




2.5 



3.5 4 

a 



4.5 



Fig. 4. Interference at the origin in a square lattice (solid curve) and lower bound (19) and upper bound (21) (dashed). 



For a = 3, 



Hz) 



4(1 + 2- 3 / 2 ) + 9(1 + 2- 5 / 2 ) + ^± + ^ I II : 



20\/2 32 \ „ , ,, 

1 WzV 

27 277 



and for a = 4, 



12 



6.5 



3.5 



0.5 



-2.5 



-5.5 




O 
X 



3.5 



6.5 



Fig. 5. Balanced TDMA scheme for square lattice for m = 3. The small dots are the lattice points Z , the crosses the transmitters, and 
the circles their receivers. Arrows indicate transmissions. The dotted lines demarcate the 3 x 3 boxes, in which there is one transmitter in 
each time slot. The transmit-receive pattern can be shifted and rotated, so that in 4 • 3 2 time slots, each node gets an opportunity to transmit 
to all four nearest neighbors with the same interference. 



Next, consider the case where the receiver is displaced diagonally, i.e., z = (r, r)/v2. In this case, we have 
4(r) = (V2 - r)~ a + (V2 + r)~ Q + 2(1 + r(r - V2))~ a / 2 + 2(1 + r(V2 + r)y a / 2 + 2(2 + r 2 )~ a / 2 
= 4(1 + 2~ a l 2 ) + a 2 (l + 2- a ' 2 - 1 ) r 2 + 0(r 4 ) . 



Since cs = c' 8 , it turns out that the excess interference coefficient is the same as for the case of axial displacement. 
D. Application to TDMA 

1) Separating the transmitters: If the transmitting nodes are located at (mZ) 2 , corresponding to an m 2 -phase 
TDMA scheme, the interference at the receiver (0,1) is I m = m~ a I(l/m), as in the one-dimensional case. 
For m > 2, the nearest 8 interferers are at the following distances: One at m — 1, two at \/T+ m 2 , two at 
y/(m — l) 2 + to 2 , one at to + 1, and two at y/Jm + l) 2 + to 2 . If these 8 terms are summed up directly and 
the Voronoi upper bound in Cor. 2 is used to upper bound the interference from the other nodes, the radius 
rb = 3to/2 — 1, since the 9-th nearest node is at distance 2m, — 1 and the Voronoi cell is to/2 wide. 

2) Balanced TDMA: As in the one-dimensional case, there is a smarter way to schedule the transmissions, such 
that the two nearest interferers are at the same distance. For all to > 1, there exists a scheme such that the eight 
nearest interferers are: Four at distance \Jm 2 + 1, two at y/ m 2 + (to — l) 2 , and two at ^to 2 + (to — l) 2 . The 
smallest distance from the receiver to a Voronoi cell of an interferer not included in these eight is 3to/2. 

We believe this balanced scheme is optimum in the sense that for all sub-lattices in which one node per to x to 
block transmits in each slot, this one causes the smallest interference. Fr om th e triangular lattice it is known that 

we cannot do better than having the six nearest interferers at distance J 2/ \/3m, and the average distance of the 
nearest six neighbors lies within 3%-5% of this bound for all to < 30. Fig. 6 illustrates the interference gain of the 
balanced scheme when compared with the simple one. 

V. Triangular lattices 

A triangular lattice, where nodes are arranged in a lattice with generator matrix G tr ; given in (2), offers the 
densest packing given the nearest-neighbor distance. It is the preferred deployment for sensor networks with isotropic 
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m 



Fig. 6. Interference for TDMA scheme for square lattices with a = 4, multiplied by m 4 for normalization. The "simple" scheme is the one 
where nodes (mis) 2 transmit, the "balanced" one is the one illustrated in Fig. 5. The dashed curves are the bounds obtained by the radial 
bounds in Cor. 2. As m — > oo, the curves converge to the same value. 



sensors, since the smallest number of sensors is needed to cover an area. The density is 2/y/3. More importantly 
though, a triangular lattice may be a good model for a CSMA-type network. Assume a high-density network, with 
A> 1 nodes per unit area, and a CSMA scheme with carrier sensing radius 1. With ideal sensing, the transmitting 
nodes cannot be denser than a triangular grid. Hence the interference in a triangular lattice is an upper bound to 
the interference in a CSMA network if transmitters are spaced at least at unit distance. 

A. Lower interference bound 

Again we can partition the interferers into rings of increasing radii; in this case, the rings are hexagons. The 
average distance to the interferers in the /c-th ring is fc(l/2+V3/4) « 0.933&, which, using Theorem 3, immediately 
yields the lower bound 



oo 



1(6) > 6 + 6 > k 



k=2 



-a 



6 + 8 'iws, 1 K*"- 1 )- 1 ' 



Lower bounding ((a) using Lemma 1, we obtain 



/ 4 \ a 2 • 3 a + 3 • 2° + 6 

/(o)>6+ v2 T7lJ 6°-i-3«-i- 2 °-i-r (22) 

B. Upper interference bounds 

Here we apply Cor. 2 with rj, = 2/ y/S, which is the smallest distance to a Voronoi cell of a point at distance y/3 
from the origin. This implies that the six nearest neighbors at distance 1 are considered separately, and all other 
points are included in the approximation. The area of a Voronoi cell is V = y/S/2, hence the bound is 

/W< S+ «. (23) 
V3 a -2 

A tighter bound is obtained by taking the nearest 18 nodes (two rings) and integrating over the Voronoi cells of 
the other points. There are six nodes at distance y/3 and six at distance 2, and r\, = y/4 + 1/3. This situation is 
illustrated in Fig. 7, where the circle has radius r&. It follows that 

J(o) < 6(1 + 2"« + 3-/ 2 ) + ^ (13/3)1 r /2 • (24) 

V 3 a — 2 

The lower bound (22) and the tighter upper bound (24) are shown in Fig. 8. 
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Fig. 7. Voronoi cells for triangular lattice. The circle of radius ^13/3 indicates the radial bound for the integration for the bound (24). 
The interference from the 18 nearest points, marked by a cross, is summed up directly, while the interference from the other nodes is upper 
bounded by integrating ||x|| _c * over the outside of the circle and multiplying by the lattice density 1/V = 2/\/3. 



20 




Fig. 8. Interference at the origin in a triangular network (solid) with lower bound (22) and upper bound (24) (dashed). The two bounds 
are uniformly tight. 



C. Transmitter-receiver offset 

Here we approximate the interference I(z) for \\z\\ < 1. As for the square lattice, we focus first on the case 
where z lies on the line between the origin and the node at (1,0), i.e., z = (||z||,0), and we let r = \\z\\. The 
interference from the 6 nearest interferers is given by 

I 6 (r) = (1 - r)" Q + (1 + r)~ a + 2(1 - r(l - r))~ a ' 2 + 2(1 + r(l + r))- Q / 2 

since the distance x of point (r, 0) to interferer (1/2, \/3/2) is given by x 2 = 1 — r(l — r). A Taylor expansion 
gives 

J 6 ( r ) = 6 + ^a 2 r 2 + 0(r 4 ) . 
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Fig. 9. Interference at position (r, 0) in a triangular network for a = 3, 3.5, 4, 5. The solid curves are the exact numerical results, while 
the dashed ones are the bounds (25). 



Taking a lower bound on the interference from the other nodes yields a lower bound for the total interference: 



I(r) > l(o) + -^a 2 r 2 



(25) 



where 7(o) denotes the lower bond on the right-hand side of (22). A better approximation may be obtained by 
decreasing the radius rj, used for the upper bounds by r, as stated in Cor. 2, but the effect of a small displacement 
r on this term is relatively small. This is illustrated in Fig. 9, where the actual interference and the quadratic lower 
bound (25) are shown. The curvature of the interference I(r) near r = 0, i.e., the excess interference, is reproduced 
quite accurately by the bound. It can further be observed that while higher path loss exponents result in smaller 
interference for no or very small displacements, there is a cross-over point after which a smaller a is better. 

Changing the direction of the displacement z does not affect the quadratic approximation, i.e., the difference 
only manifests itself for r > 1/4. So in the triangular case, the offset coefficient c cx ~ 3a 2 /2 for all directions. 



D. Application to TDMA 

1) Separating the transmitters: The most straightforward TDMA scheduling scheme for the triangular lattice 
is probably the one where one node in a rhombus containing m 2 nodes transmits. The generator matrix for the 
transmitting lattice is then 

"l{modd}/2 
mv/3/2 



G 



rh 



m 




For even m, the transmitters form a rectangular lattice with horizontal spacing m and vertical spacing m\/3/2. If 
the receiver is to the right of each transmitter, the interference is acceptable, but this is certainly not the optimum 
scheme, since the nearest interferer is at distance m — 1. An improved scheduler lets one node in a parallelogram 
containing (m + l)m nodes transmit: 

m + 1 —m/2 
mV3/2 



G 



par 



2 ) Balanced TDMA: An even better scheme is the one inspired by the balanced one-dimensional scheme, where 
in each horizontal row of nodes, the pattern ( T^R m ~ 1 ) m is employed. This way, the four closest interferers are all 
at distance m, and the next two at y/ (m/2 + l) 2 + 3m 2 /4. Since these distances are close to the sphere-packing 
bound, this may be the optimum scheme. 

The three TDMA MACs are compared in Fig. 10. 
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Fig. 10. Interference for TDMA schedulers in the triangular lattice for a = 4. For a fair comparison, the area of the corresponding Voronoi 
cells is used as the x axis, and the y axis is scaled by \~ a/2 = A~ 2 , where A = l/V, since scaling the network by a factor s in both 
dimensions scales the density by s~ 2 and the interference by s~ a . 



VI. Conclusions 

By judicious application of Jensen's inequality, tight and general upper and lower bounds on the interference in 
lattice networks can be calculated. The lower bounding technique yields expressions involving the Riemann zeta 
function, which, in turn, can be tightly lower bounded in closed-form (Lemma 1). 

When transmitters are arranged in a lattice, but their receivers are not, such as in a cellular downlink scenario 
or when a MAC scheme focuses on separating the transmitters, the interference needs to be characterized at some 
distance r from the desired transmitter. We have introduced a simple yet effective quadratic approximation that 
shows how the interference increases with growing link distance r. In terms of the signal-to-interference-and-noise 
ratio, increasing r has two negative effects: it reduces the strength of the desired signal, and it increases the 
interference. Our framework captures both. 

For the case where all nodes in a network form a one-dimensional, square, or triangular network, we have 
analyzed the interference induced by basic TDMA schemes, and we have suggested superior balanced schemes 
that are optimum or near-optimum. It turns out that scheduling transmitter-receiver pairs instead of just transmitters 
significantly enhances the throughput. 

The interference bounds are also useful to avoid long simulation times. Normally as a \. 2 (in two-dimensional 
networks), the diameter of the simulated networks has to be impractically large just to get a good approximation 
of the total interference. The upper bounds obtained from integration over Voronoi cells (Theorem 1) demonstrate 
that it is possible, with negligible loss in accuracy, to assume that the interference from nodes outside a certain 
radius is the same as the expected interference in a network where these nodes are distributed as a homogeneous 
Poisson point process. 
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